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Abstract 



IS 



A measurement of the Z/j* transverse momentum (p^) distribution in proton-proton colhsions at ^/s ~ 7 TeV 
presented using Z/j* — > e+e^ and Z/j* — > /i+/^~ decays collected with the ATLAS detector in data sets with integrated 
luminosities of 35 pb~^ and 40 pb~^, respectively. The normalized differential cross sections are measured separately 
for electron and muon decay channels as well as for their combination up to p^ of 350 GeV for invariant dilepton 
masses 66 GeV < ma < 116 GeV. The measurement is compared to predictions of perturbative QCD and various event 
generators. The prediction of resummed QCD combined with fixed order perturbative QCD is found to be in good 
agreement with the data. 
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In hadron collisions, the weak vector bosons W and Z 
are produced with a momentum component transverse to 
the beam axis, which is balanced by a recoiling hadronic 
system mainly arising from initial state QCD radiation of 
quarks and gluons. The measurement of the boson trans- 
verse momentum offers a very sensitive way of studying dy- 
namical effects of the strong interaction, complementary to 
measurements of the associated production of the bosons 
with jets [1]. The simple signatures of Z/j* — ^ e+e~ 
and Z/j* — )• /i^/i^ production, which can be identified 
with little background, enable a precise measurement of 
the boson transverse momentum (pf ) and thus provide 
an ideal testing ground for predictions of QCD and phe- 
nomenological models. Moreover, the knowledge of the p^ 
distribution is crucial to improve the modelling of W bo- 
son production needed for a precise measurement of the 
W mass [2, 3], in particular in the low p^ region which 
dominates the cross section. 

This Letter presents a measurement of the Z/^* 
normalized transverse momentum distribution in proton- 
proton collisions at y/s — 7 TeV using Z/j* — >• e+e~ 
and Z/'y* — >■ jJ.^ ^^ decays collected in 2010 with the 
ATLAS detector in data sets with integrated luminosities 
of 35 pb~^ and 40 pb~^, respectively. In the normalized 
transverse momentum distribution, many systematic 
uncertainties cancel. In particular, the precision of the 
measurement is not impaired by the uncertainty on 
the integrated luminosity. The normalized transverse 
momentum distribution is defined in the following way: 
1/(7^'^ X dcr'^'^/dpf, where a^'^ is the measured inclu- 
sive cross section of pp — > Z/"f* + X multiplied by the 
branching ratio of Z/^* -> £+£^ within the fiducial 
acceptance, and X denotes the underlying event and 



the recoil system. For both the ee and p,p, channels, the 
fiducial acceptance is defined by the lepton transverse 
momentum and pseudorapidity^ , and by the invariant 
mass of the lepton pair m«: p^ > 20 GeV, \rj^\ < 2.4 
and 66 GeV < to« < 116 GeV. The measurements in 
both decay channels are corrected for detector effects and 
QED final state radiation (FSR). The combined result is 
compared to predictions of perturbative QCD calculations 
and QCD-inspired models implemented in various event 
generators. 

Predictions and previous measurements are discussed 
in Section 2. The ATLAS detector and trigger are de- 
scribed in Section 3. In Sections 4 and 5 the event simula- 
tion and selections are described. Section 6 reports the p^ 
measurements for different treatments of QED final state 
radiation. Systematic uncertainties are discussed in Sec- 
tion 7. Section 8 presents the combined result which is 
compared with various models. 

2. QCD predictions and previous measurements 

Perturbative QCD (pQCD) calculations have been per- 
formed up to ©(ag) in the strong coupling constant as 
[4, 5] and are expected to be reliable at large pf . In this 
kinematic regime, the cross section is dominated by the ra- 
diation of a single parton. Fully differential inclusive boson 
production cross sections can be obtained at ©(ag) with 



^The nominal pp interaction point at tiie centre of the detector 
is defined as the origin of a right-handed coordinate system. The 
positive a;-axis is defined by the direction from the interaction point 
to the centre of the LHC ring, with the positive y-axis pointing up- 
wards. The azimuthal angle ip is measured around the beam axis and 
the polar angle 9 is the angle from the z-axis. The pseudorapidity is 
defined as 77 = — lntan(6/2). 
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Fewz [6, 7] and Dynnlo [8]. While the integrated 0(a§) 
cross section predictions are finite, the fixed order pQCD 
prediction diverges at vanishing p^. In this regime, the 
leading contribution of multiple soft gluon emissions to the 
inclusive cross section can be resummed to all orders [9- 
11] up to next-to-next-to-leading logarithms (NNLL) [12] 
in as- The Resbos generator [13] matches the prediction 
of soft gluon resummation including a non-perturbative 
form factor [14] at low pf with the fixed order pQCD cal- 
culation at 0{as) at high pf , which is corrected to 0{ag) 
using ivT-f actors. 

Similar to resummed calculations, parton showers 
provide an all-order approximation of parton radiation 
in the soft and collinear region. In order to describe the 
large pf region, the parton shower based leading-order 
event generators Pythia [15] and Herwig [16] apply 
weights to the first or hardest branching, respectively, to 
effectively merge the 0{ag) and 0{as) pQCD predictions. 
The next-to-leading order (NLO) Monte Carlo generators 
Mc@NLO [17] and Powheg [18] incorporate NLO QCD 
matrix elements consistently into the parton shower 
frameworks of Herwig or Pythia. 

The Alpgen [19] and Sherpa [20] event generators 
implement tree-level matrix elements for the generation of 
multiple hard partons in association with the weak boson. 
The matrix-element calculations for various parton mul- 
tiplicities are matched with parton showers (which in the 
case of Alpgen are provided by either Pythia or Her- 
wig) such that double counting is explicitly avoided by 
means of weighting procedures [21] or veto algorithms [19]. 

As the predictions of these generators differ signifi- 
cantly and show a considerable dependence on adjustable 
internal parameters [22] , a precise measurement of the bo- 
son transverse momentum distribution is an important in- 
put to validate and tune these models. 

The Z/"f* boson px distribution has been measured 
in proton-antiproton collisions at the Tevatron collider at 
centre of mass energies of y/s = 1.8 TeV and 1.96 TeV 
[23-27]. For p^ < 30 GeV, these measurements found a 
good agreement with Resbos and disfavoured models [28] 
suggesting a broadening of the p^ distribution for small x 
values [25, 27], where x is the fraction of the momentum 
of one of the two partons with respect to the proton mo- 
mentum. At large p^, the ©(ckg) pQCD prediction was 
reported to underestimate the measured cross section by 
up to about 25% [25, 26]. 

3. The ATLAS detector 

The ATLAS detector system [29] comprises an inner 
tracking detector immersed in a 2 T axial magnetic field, 
a calorimeter, and a large muon spectrometer with a su- 
perconducting toroid magnet system. Charged particle 
tracks and vertices are reconstructed with silicon pixel 
and strip detectors covering \r]\ < 2.5 and transition ra- 
diation detectors covering \r]\ < 2.0. These tracking de- 
tectors are surrounded by a finely segmented calorimeter 



system which provides three-dimensional reconstruction of 
particle showers up to j?]] < 4.9. The electromagnetic com- 
partment uses liquid argon as the active material and is 
divided into barrel {\r]\ < 1.5), end-cap (1.4 < [jy] < 3.2) 
and forward (3.1 < \r]\ < 4.9) components. The hadron 
calorimeter is based on scintillating tiles in the central re- 
gion {\r]\ < 1.7). It is extended up to 4.9 in pseudorapidity 
by end-caps and forward calorimeters which use liquid ar- 
gon. The muon spectrometer is based on three large super- 
conducting toroids arranged with an eight-fold azimuthal 
coil symmetry around the calorimeters, covering a range of 
\r]\ < 2.7 and providing an integral magnetic field varying 
from about 1 to 8 T • m. Three stations of drift tubes and 
cathode strip chambers enable precise muon track mea- 
surements, and resistive-plate and thin- gap chambers pro- 
vide muon triggering capability and additional measure- 
ments of the if coordinate. 

The ATLAS detector has a three-level trigger system 
which reduces the event rate to approximately 200 Hz be- 
fore data transfer to mass storage. The triggers employed 
require the presence of a single electron or muon candidate 
with pt > 15 GeV or px > 13 GeV, respectively. Lower 
thresholds were used for the early data. The trigger effi- 
ciencies are defined for Z/7* — >■ e+e" and Z/j* -^ f^'^IJ-' 
events as the fraction of triggered electrons or muons with 
respect to the reconstructed lepton and are studied as a 
function of their px and 77. For muons, the efficiencies 
are also obtained separately in ip regions which match the 
geometry of the trigger chambers. The efficiency for sin- 
gle leptons is derived from data using Zf-f* -^ t^t' can- 
didate events or using independent triggers by matching 
reconstructed lepton candidates to trigger signals in the 
calorimeter (muon spectrometer) in case of the ee (^[ip) 
decay channel. For pi^ > 20 GeV, the efficiency is 99% for 
electrons and 77% (93%) for muons in the barrel (end-cap). 
For signatures with two high-Sx electrons, the trigger is 
fully efficient. The trigger efficiency for Z/j* — ?> /^"""M" 
events is determined to be on average 97.7% and to be con- 
stant as a function of p^ within an uncertainty of 0.1— 0.7% 
depending on p^ . 

4. Event simulation 

The properties, including signal efficiencies and 
acceptances, of Z/'y* — >■ e+e~, Z/'y* — >■ n'^fj,^ and back- 
ground processes are modelled with Pythia [15] using 
the MRST2007LO* [30] parton distribution functions 
(PDF), Mc@NLO [17] and Powheg using CTEQ6.6 [31] 
PDFs. Mc@NLO uses Herwig for the parton shower 
and Jimmy [32] for the underlying event. Powheg is 
interfaced to Pythia for the underlying event and the 
parton shower. The event generators are interfaced to 
Photos [33] to simulate QED FSR. Version 6.4 of Pythia 
is used with the px-ordered parton shower and with pa- 
rameters describing the properties of the underlying event 
which were tuned to Tevatron measurements [34]. For 
systematic studies and comparisons, a Mc@nlo based 



signal sample is used with underlying event parameters 
(Jimmy) tuned to Tevatron and 7 TeV ATLAS pp collision 
data [35]. The response of the ATLAS detector to the 
generated particles is modelled using Geant4 [36], and 
the fully simulated events [37] are passed through the 
same reconstruction chain as the data. The Monte Carlo 
simulation (MC) is corrected for differences with respect 
to the data in the lepton reconstruction and identifica- 
tion efficiencies as well as in energy (momentum) scale 
and resolution. The efficiencies are determined from a 
tag-and-probe method based on reconstructed Z and W 
events [38, 39], while the resolution and scale corrections 
are obtained from a fit to the observed Z boson line 
shape. The lepton identification efficiencies can depend 
on the hadronic activity, which is correlated with the 
Z/"f* transverse momentum. Therefore, using the tag- 
and-probe method, it is verified that the p^ dependence 
of the single lepton efficiency is correctly modelled after 
efficiency corrections. Differences between data and sim- 
ulation are mostly consistent with statistical fluctuations 
and are considered as systematic uncertainties due to the 
modelling of the efficiencies as described in Section 7. 

Multiple pp interactions per bunch crossing (pileup) 
are accounted for by overlaying simulated minimum bias 
events. To match the observed instantaneous luminosity 
profile, the MC events are reweighted to yield the same 
distribution of the number of primary vertices as measured 
in the data. 

In the following, if not stated otherwise, the gener- 
ated samples are fully simulated, pileup and resolution 
corrected, and interfaced to Photos, with Pythia as 
the default MC signal sample. The background samples 
are produced with Pythia for W — >■ £v, with Mc@NLO 
and POWHEG for tt, and with Pythia and Mc@nlo for 

Z/7* -^T+T-. 

5. Event selection 

The analysis uses data taken during stable beam con- 
ditions with properly operating inner detector, magnets, 
and calorimeter or muon spectrometer, in case of the ee or 
/u/z channel, respectively. Events are required to have at 
least one primary vertex reconstructed from at least three 
tracks. Z/Y events are selected by requiring two oppo- 
sitely charged electrons or muons, defined below, with an 
invariant mass 66 GeV < mu < 116 GeV. 

Electrons are reconstructed from the energy deposits 
in the calorimeter matched to inner detector tracks. They 
are required to have a transverse energy E^ > 20 GeV and 
pseudorapidity \ri'^\ < 2.4, excluding the transition regions 
between the barrel and end-cap calorimeter components at 
1.37 < \r)'^\ < 1.52. They should pass medium identifica- 
tion criteria based on shower shape and track quality vari- 
ables [38] to provide rejection against hadrons. The single 
electron selection efficiency varies in the range 90 — 96% 
depending on pseudorapidity and azimuthal angle. 



Muons are reconstructed from matching tracks in the 
inner detector and muon spectrometer with pt^ > 20 GeV 
and [77'' I < 2.4 measured using the information from these 
two detector subcomponents. To ensure the compatibility 
of the muon track with the primary vertex, the correspond- 
ing impact parameters in the transverse and longitudinal 
direction with respect to the beam axis have to be smaller 
than 1 mm and 5 mm, respectively. The muon candi- 
dates are required to be isolated to suppress background 
from heavy flavour production using the transverse mo- 
mentum sum of tracks with p^ > 1 GeV within a cone 
of size AR = ^/{Ar])'^ + {Atp)"^ = 0.2 around the muon 
track. This sum has to be smaller than 0.2 x pi^. The 
single muon selection efficiency is of the order of 87 — 95% 
depending on pseudorapidity and azimuthal angle. 

After this selection, 8 923 Z/j* -^ e+e" and 15 060 
Z/j* — > i-i'^ IJ-'^ candidate events are found in the data. 
The main reasons for the difference in the number of can- 
didate events in the ee and fifi decay channels are: the 
integrated luminosity for which all relevant detector com- 
ponents were operating properly; regions in the electro- 
magnetic calorimeter with readout problems; a reduced 
acceptance for electrons in the transition region between 
the barrel and end-cap calorimeter. Figure 1 demonstrates 
the agreement of the data and the simulation in the dilep- 
ton mass spectrum for the selected events including the 
background, after applying the corrections described in 
Section 4. 

The background contribution from Z/"f* — > t^t^, W 
boson, and tt production is estimated as a function of p^ 
using simulation, where the cross sections are normalized 
to next-to-next-to- leading order (NNLO) predictions for 
Z/7* and W and NLL-NLO predictions for tt production. 
The procedure outlined in Ref. [38] is followed here. 

For both ee and fifi channels, the main background at 
high pf arises from tt production. At low p^ the back- 
ground is dominated by QCD multijet production, where 
a jet is falsely identified as a primary e or fi. Its contribu- 
tion is determined from the data as follows. 

In the ee channel, the normalization of the multijet 
contribution is derived from a fit of signal and background 
templates to the observed dilepton invariant mass distri- 
bution with loosened identiflcation requirements for one 
of the two reconstructed electron candidates. An extended 
mass range of 50 GeV < niee < 130 GeV is used which pro- 
vides a better background constraint in the off-resonance 
region. The normalization derived from this loosened se- 
lection has to be scaled to the Z/'f* event selection, which 
requires two reconstructed electron candidates of medium 
quality. The scaling factor is determined from a QCD 
multijet enhanced control sample with single electron can- 
didates which fulfil the loosened electron identification re- 
quirements. The contamination with other events, in par- 
ticular the contribution from W production, is suppressed 
by rejecting events with large missing transverse energy. 
The remaining contamination is determined using simu- 
lated events. The systematic uncertainty of the normal- 
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Figure 1: The observed (a) dielectron and (b) dimuon invariant mass distributions compared to simulation. The shaded bands indicate the 
systematic uncertainty on the MC prediction due to the lepton momentum resolution correction. The Pythia sample, which is normalized 
to the next-to-next-to-leading order (NNLO) prediction [40, 6, 7], is used for signal events. The total background, which is invisible on this 
scale, amounts to only 1.5% and 0.4% in the ee and fifi decay channels, respectively. 



ization is determined by varying the background templates 
and the criteria for the loosened selection. The shape 
of the multijet background as a function of p^ is deter- 
mined from a dielectron sample with an invariant mass 
66 GeV < rriee < 116 GeV for which exactly one elec- 
tron passes and one fails the medium identification criteria. 
The difference between same and opposite sign events is 
taken as the shape uncertainty. 

For the fj,fi selection, the multijet contribution is es- 
timated from four two-dimensional regions which are ob- 
tained by changing the mass window to 40 GeV < m^^ < 
60 GeV or by inverting the isolation criterion. The nor- 
malization is determined from the number of candidate 
events in these four regions which is corrected by the ex- 
pected number of non-QCD multijet background events 
and number of signal events. The number of signal events 
is determined from the measurement in the signal region 
which is corrected for the background including the mul- 
tijet background, and extrapolated to the other regions 
using relative efficiencies extracted from the simulation. 
The resulting equation can be solved for the number of 
QCD multijet background events. The shape is deter- 
mined from the control region where both muon candi- 
dates are non-isolated and their invariant mass is within 
66 GeV < TTi^p < 116 GeV. Systematic uncertainties are 
determined by comparing the results to those of alterna- 
tive methods which use template fits or same sign and 
opposite sign events. 

The total background contribution, which amounts to 
(1.5 ± 0.6)% and (0.4 ± 0.2)% in the ee and fifi chan- 



nels, respectively, is found to increase from 0.5% (0.2%) 
to 3.5% (1.5%) as a function of p^ in the ee (nfJ.) analysis. 

6. The measured p^ distribution 

The Z/"f* transverse momentum is reconstructed from 
the measured lepton momenta. The p^ range is divided 
into 19 bins from GeV to 350 GeV with widths of 3 GeV 
for pf < 30 GeV and increasing widths at larger trans- 
verse momenta as given in Table 1. For this binning, the 
fraction of simulated events reconstructed in a particular 
p^ bin which have generator-level p^ in the same bin is 
always better than 60% and reaches values above 90% in 
the highest pf bins for both decay channels. The bin- 
by-bin efficiency, defined as the ratio between the number 
of signal events which pass the final selection and the to- 
tal number of generated events within the fiducial region, 
are on average about 56% in the ee and 83% in the /z^ 
channel, respectively. Three Z/j* candidate events with 
pf > 350 GeV are found, which are not considered further 
due to the limited statistical significance. 

The observed pf spectrum is found to be well described 
by the simulation, using the default MC signal samples 
and background estimations as described in Sections 4 
and 5. A bin- by-bin efficiency correction is used to correct 
(unfold) the observed data for detector effects and QED 
FSR, where the correction factors are determined from the 
default MC signal sample. Alternative matrix unfolding 
methods [41, 42], which explicitly take the bin-to-bin mi- 
gration into account, yield compatible results. However, 



these techniques require higher data statistics to fully ex- 
ploit their advantages. Therefore they are currently used 
only to estimate the systematic uncertainties due to the 
unfolding method as discussed in Section 7. 

In Table 1, the cross section measurements in the ee 
and /i/^ decay channels are reported in the fiducial volume, 
which is defined by the lepton acceptance p^ > 20 GeV 
and \ri^\ < 2.4, and the invariant mass of the lepton pair 
66 GeV < m« < 116 GeV. This imphes for the ee de- 
cay channel a small acceptance correction due to the dis- 
carded events, in which one or more electrons are within 
the calorimeter transition region, 1.37 < \ri^\ < 1.52. The 
resulting correction of the normalized differential cross sec- 
tion is smaller than 0.6%. The measurement is reported 
with respect to three distinct reference points at particle 
level regarding QED FSR corrections. The true dilepton 
mass mu and transverse momentum p^ are either defined 
by the final state leptons after QED FSR ("bare" leptons), 
or by recombining them with radiated photons within a 
cone of AR = 0.1 ("dressed" leptons), or by the Z/j* 
propagator. The propagator definition corresponds to a 
full correction for QED FSR effects, allows for a combi- 
nation of the electron and muon channels, and facilitates 
a direct comparison of the measurement with QCD calcu- 
lations. The QED FSR corrections are at most 8% (5%) 
for the normalized differential cross section in the ee (/i/i) 
decay channel. 

7. Systematic uncertainties 

Systematic uncertainties arise mainly from lepton 
efficiencies, momentum scale and resolution, and from 
the unfolding procedure. They are evaluated by varying 
separately each parameter in question and recalculating 
the bin-by-bin correction factors used to unfold the 
p^ spectrum. The observed deviation in the unfolded 
distributions is limited by the statistics of the simulated 
events. The statistical component of these deviations 
is estimated using a bootstrap method [44] based on 
multiple resampling of the simulated events. Each sample 
has the same size and may contain the same events 
multiple times. Symmetric systematic uncertainties are 
derived in such a way [45] that the area spanned by these 
uncertainties covers 68% of the integral over a Gaussian 
with mean and width equal to the mean and standard 
deviation of the bootstrap distributions. The resulting 
uncertainties are equal to the standard deviation if the 
bootstrap distribution is centred at zero, and approach 
the mean plus half the standard deviation if the mean 
deviates from zero. 

The following sources of systematic uncertainties on 
the measured normalized differential cross section are eval- 
uated, where the quoted uncertainties are relative percent- 
ages. 

Lepton reconstruction and identification: The efficien- 
cies are determined as a function of p^ using a tag-and- 
probe method. Due to the limited statistics at high p^, 



the uncertainties are parametrized by a function which 
increases with p^. For the ee analysis, the uncertainty 
due to the electron reconstruction and identification varies 
between 1.0% and 3.2%. An additional uncertainty of 
0.1% arises from the modelling of local calorimeter readout 
problems. In the case of the /i/x selection, the uncertain- 
ties due to the efficiencies on the trigger, on the muon 
reconstruction and identification, and on the isolation re- 
quirement are evaluated separately and are found to be 
within 0.4%, 2.3%, and 1.0%, respectively, except for the 
three highest pf bins, where they reach 0.8%, 4.9%, and 
2.7%. 

Lepton energy (momentum) scale and resolution: The 
scale and resolution corrections of the simulation are var- 
ied within their uncertainties estimated from the fit to the 
observed Z/7* line shape. Correlations across pf bins are 
taken into account. Due to the normalization to the inclu- 
sive cross section, a systematic shift at low p^ is balanced 
by a shift in the opposite direction at high p^. In the 
ee {fj.fl) analysis, the uncertainty due to scale variations 
is found to be 2.7% (0.2%) for the lowest bin, decreas- 
ing down to 0.2% (< 0.1%) at p| - 10 GeV, and then 
increasing up to 4.4% (0.4%) at the highest p^ values. 
Uncertainties due to resolution are estimated to be 0.5% 
for the ee channel and between 0.1% and 0.7% for the fifi 
channel. 

Unfolding procedure: The bin-by-bin correction factors 
depend on the shape of the assumed underlying p^ dis- 
tribution, which leads to a systematic uncertainty evalu- 
ated in the following way. The default MC signal sample 
(Pythia) is reweighted to different true p^ shapes using 
Resbos and Mc@nlo. The variance of these generator 
predictions does not entirely cover the observed differ- 
ence between simulation and data. Therefore the Pythia 
signal sample is reweighted, in addition, to distributions 
based on unfolded data. The spectra obtained by unfold- 
ing the data either with the bin-by-bin method or alter- 
native matrix unfolding techniques [41, 42] are considered. 
These new corrected spectra feature different uncertain- 
ties. The bin-by-bin method suffers a larger systematic 
uncertainty due to the assumed true p^ shape, whereas 
matrix-based unfolding is nearly independent of this as- 
sumption, but suffers from a larger statistical uncertainty. 
Each of the reweighted spectra is treated in the same way 
as the data, and is unfolded with the default bin-by-bin 
correction factors. The maximum deviation from the re- 
spective true p^ spectrum is considered as a systematic 
uncertainty. A possible infiuence from the parton shower 
and hadronization model is estimated by comparing the 
bin-by-bin correction factors determined from either the 
default Pythia sample or the Mc@nlo sample. Mc@nlo 
uses Herwig for the parton shower and Jimmy for the 
underlying event. To separate these model uncertainties 
from the uncertainty due to the underlying p^ distribu- 
tion, which is already accounted for, the Mc@nlo sample 
is reweighted to the pf shape of the default MC signal sam- 
ple. The uncertainties due to the shape of the p^ distri- 



Table 1; The measured normalized differential cross section 1/a^'^da^'^/dp^ in bins of p^ for Z/'y* — > e+e~ and Zf-y* — > ^i~^li~ production. 
The cross sections, which are to be multiplied by the factor k, are reported with respect to three different treatments of QED final state 
radiation for the definition of lepton and Z/'y* momentum at particle level: corresponding to the Z/'y* propagator (propag.), to dressed 
leptons (with recombined radiated photons within a cone of AR = 0.1), and to bare Icptons, respectively. The relative statistical (stat.) and 
total systematic (syst.) uncertainties are given. The data can be obtained electronically through the HepData repository [43]. 
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10-2 
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7.22 


7.24 


7.26 
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6.45 


6.46 


6.45 
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4.89 


4.88 


4.85 
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4.63 
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3.66 


3.64 


3.59 
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3.94 
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3.26 


3.25 


3.20 
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2.88 


2.86 
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6.2 


2.6 
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2.14 


2.13 


2.08 


10-3 


9.2 


5.4 


2.14 


2.13 


2.11 


10-3 


7.2 


2.7 
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1.21 


1.20 


1.17 


10-3 
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1.31 


1.30 


1.28 
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3.0 
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5.69 


5.63 


5.44 
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5.52 


5.47 


5.40 


10-4 
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1.74 


1.73 


1.67 


10-4 


9.6 


6.1 


1.52 


1.51 


1.49 


10-4 


7.5 


4.4 


180-350 


0.78 


0.77 


0.73 


10-5 


27.0 


7.8 


1.14 


1.14 


1.11 


10-5 


18.9 


6.6 



bution and due to the modelling of the parton shower and 
the hadronization are combined to yield the total unfold- 
ing uncertainty, which is found to be within 2.0% (1.3%) 
in the ee (/x/x) channel for p^ between 6 GeV and 100 GeV. 
For pf < 6 GeV the uncertainty is as large as 3.6% (4.7%) 
and for pf > 100 GeV it is as large as 4.2% (2.9%) in the 
ee (fJ-fJ.) channel. The unfolding uncertainty is dominated 
by the deviations observed when reweighting the default 
MC signal sample to the pf distributions obtained from 
the data with the matrix unfolding techniques. 

Background contamination: Uncertainties in the esti- 
mation of the background from QCD multijct, weak boson, 
and tt production yield values of up to 1.4% (0.6%) for the 
ee (fifi) analysis when propagated to the normalized dif- 
ferential cross section. 

Modelling of pileup corrections: Pileup has a small in- 
fluence on this measurement. An uncertainty of 0.3% on 
the normalized differential cross section is derived. 

MC sample statistics: The uncertainties are within 
0.4% - 1.5% (0.3% - 0.8%) in the ee {pp) channel, ex- 
cept for highest p^ values for which the uncertainties reach 
3.6% (1.6%). 

QED final state radiation: A conservative systematic 
uncertainty of 0.6% due to the p^-dependent modelling of 
QED FSR is assigned. This uncertainty addresses both 



potential differences between the approximation used in 
Photos compared to exact second order QED FSR ma- 
trix element calculations [46, 47] and uncertainties in the 
simulation of the interaction of the radiated photons with 
the detector material [38]. 

The systematic uncertainties listed above are added 
quadratically to obtain the total systematic uncertainties 
listed in Table 1. 

8. Results and conclusions 

The normalized differential cross section measurements 
for Z/^* — > e^e- and Z/j* — > p'^fJ.^ production are in 
good agreement with each other at the Z/j* propagator 
level; see Figure 2 for the ratio of the measured cross sec- 
tions. The two decay channels are combined at Z/^* prop- 
agator level using a x^ minimization method which takes 
into account the correlated systematic uncertainties for the 
ee and pp channels [48]. The uncertainties due to the 
unfolding procedure and QED FSR are considered to be 
common for the two channels. The minimization yields a 
x2/d.o.f. = 17.0/19 indicating the excellent compatibility 
of the electron and muon data. 

The combined measurement of the normalized differen- 
tial cross section within the fiducial lepton acceptance as a 
function of pf , 1/cr^'^da^^ /dp^, is shown in Figure 3 and 
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Figure 2: 

Z/7* — >• e ' e ana z,/y — s> /^ ' /i 

Z/7* propagator level using an integrated luminosity of 35 pb~^ and 

40 pb~^, respectively. The error bars shown include statistical and 

systematic uncertainties. The systematic uncertainties due to the 

unfolding procedure and QED FSR, which are correlated between 

the electron and muon decay channel, arc omitted. 



Table 2. In addition, the acceptance corrections Ac needed 
to extrapolate the measurement to full lepton acceptance, 
but keeping the mass range 66 GeV < m« < 116 GeV, 
are reported. They are determined using Pythia and the 
MRST2007LO* PDF set. The acceptance for the inclusive 
fiducial cross section is 0.48. However, the acceptance cor- 
rections Ac for the normalized differential cross section are 
within 10% of 1.0 for the bins with pf < 80 GeV. The un- 
certainty on Ac is estimated by: reweighting the Pythia 
prediction using the HERAPDFl.O [49] and CTEQ6.6 [31] 
parton distribution functions; propagating the CTEQ6.6 
PDF error eigenvector sets; and by taking into account 
the difference to the predictions obtained with Mc@nlo, 
Resbos, and Fewz. 

In Figure 4, the measurement within the fiducial ac- 
ceptance is compared with predictions of pQCD calcula- 
tions and of various event generators introduced above. 
The ©(as) and 0{ag) pQCD predictions of the pf depen- 
dent cross section are obtained with Fewz v2.0 [7] and 
the MSTW2008 PDF sets [50]. The inclusive cross sec- 
tion, which is used to normalize the prediction, is calcu- 
lated in the same way. The uncertainties on the normal- 
ized predictions are evaluated by variation of the rcnormal- 
ization and factorization scales by factors of two around 
the nominal scale fiR — fxp — Mz with the constraint 
0.5 < fJ-B./tJ-F ^ 2, by variation of as within a range corre- 
sponding to 90% confidence-level limits [51], and by using 
the PDF error eigenvector sets at 90% confidence- level. 



Table 2: The combined normalized differential cross section at Z/7* 
propagator level, 1/a da /dp^, as a function of the average Z/'y* 
transverse momentum (p^) with relative statistical and total system- 
atic uncertainties. The multiplication with the inverse acceptance 
correction A^ (given with uncertainties, "unc") yields the nor- 
malized differential cross section l/(T*°'do-'°'/dp^ extrapolated to 
the full lepton acceptance. The data can be obtained electronically 
through the HepData repository [43]. 
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They amount to ^ 10% and 



for the ©(as) and 



©(ag) prediction, respectively, with a dominant contribu- 
tion of 9% and 6.5% from the scale variations. In contrast 
to the Z/j* inclusive cross section, the prediction of the 
p^ distribution suffers from substantial scale uncertainties 
indicating non-negligible missing higher order corrections. 
For p^ > 18 GeV, the pQCD prediction receives an ©(ag) 
correction of 26 — 36%. Despite this correction, the ©(ag) 
prediction undershoots the data by about 10%, which is 
comparable to the size of the scale uncertainty. This deficit 
is smaller compared to the 15 — 25% difference observed at 
the Tevatron [25, 26]. At low boson transverse momenta, 
where fixed order pQCD calculations are not expected to 
give an adequate description of the cross section, the dis- 
agreement increases rapidly towards vanishing p^. 

In addition, the measurement is compared to the pre- 
dictions of Resbos and various event generators. The con- 
sistency with the data is verified with a x^ test which uses 
the x^ definition also used for the combination of the ee 
and fj,fi decay channels. 

The Resbos [13] prediction, which combines re- 
summed and fixed order pQCD calculations, is based on 
the CTEQ6.6 [31] PDF set and a resummation scale of 
mz- It is verified that the different PDF sets used for the 
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Figure 3: The combined normalized differential cross section at Z/7* propagator level as a function of p^ for (a) the range p^ < 30 GeV 
and (b) the full range compared to the predictions of Resbos, Pythia, and Fewz at 0{a'^). The error bars shown include statistical and 
systematic uncertainties. For the combination, the ee {nfJ.) channel contributes with an integrated luminosity of 35 pb~^ (40 pb~^). At low 
p^ the Fewz prediction diverges and is omitted. 
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Figure 4: Ratios of the combined data and various predictions over the Resbos prediction for the normalized differential cross section as 
a function of p^: (a) Fewz predictions at 0{as) and 0{ag); (b) predictions from the generators Pythia, McOnlo, Powheg, Alpgen and 
Sherpa. The Fewz predictions are shown with combined scale, as, and PDF uncertainties. The data points are shown with combined 
statistical and systematic uncertainty. At low p^ the 0{as) and 0{a'g) predictions of Fewz diverge and are omitted. 



Fewz and Resbos predictions lead to differences below 
3%. Resbos shows good agreement with the measurement 
over the entire pf range (x^/d.o.f. = 21.7/19), indicating 
the importance of resummation even at relatively large 
p^. However, its predictions are slightly higher than the 
data for p^ values in the range of 10 GeV to 40 GeV and 
slightly lower above 40 GeV. 

The Alpgen [19] and Sherpa [20] generators consider 
processes with up to five additional hard partons associ- 
ated with the produced boson and give a good descrip- 
tion of the entire measured spectrum, up to large pf , with 
X^/d.o.f. of 31.9/19 and 16.8/19, respectively Here, the 
enhancement of the cross section compared to the 0(ag) 
prediction can be attributed to processes with large parton 
multiplicities [52], which correspond to tree- level diagrams 
of higher order in the strong coupling. Sherpa vl.2.3 and 
Alpgen v2.13 are used, with the latter being interfaced to 
Herwig v6.510 [16] for parton shower and fragmentation 
into particles, and to Jimmy v4.31 [32] to model underly- 
ing event contributions. For Alpgen, the CTEQ6L1 [53] 
PDF set is employed and the factorization scale is set to 
H% = mjj + X]Pt' where the sum extends over all associ- 
ated partons. The Sherpa prediction uses the CTEQ6.6 



PDF set and /i|n 



(pf) 



The predictions of the parton shower event generators 
Pythia and Mc@nlo are based on the simulated samples 
as described above. Figure 4b also shows the predictions of 
Powheg vl.O [18] interfaced to a Pythia version with an 
underlying event tune to Tevatron and 7 TeV pp collision 
data [54]. Whereas Mc@nlo (xVd-o.f. = 111.6/19) and 
Powheg (x^/d.o.f. = 100.4/19) deviate from the data 
at low and high p;f , Pythia describes the measurement 
well over the entire range of boson transverse momentum 
(xVd.o.f. = 17.9/19). 

In summary, the Z /^* transverse momentum differen- 
tial distribution has been measured up to pf — 350 GeV 
for electron and muon pairs with invariant masses 
66 GeV < mu < 116 GeV produced in pp collisions at 
-y/i = 7 TeV based on integrated luminosities of 35 pb~^ 
and 40 pb"-'^, respectively, recorded with the ATLAS de- 
tector. Resbos describes the spectrum well for the entire 
pf range. At pf > 18 GeV, the central Fewz ^^(ag) pre- 
diction underestimates the data by about 10%, which is 
comparable to the size of the combined experimental and 
theoretical uncertainty. The measurement is compared 
to predictions of various event generators and a good 
agreement with Sherpa, Alpgen, and Pythia is found. 
Except for the lowest p^ values, the measurement is 
limited by statistics rather than systematic uncertainties. 
The systematic uncertainties are also mostly limited by 
the size of the data sample and are expected to improve 
with increasing integrated luminosity. 
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